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I .  Background 

The  objective  of  this  research  effort  has  been  to  develop 
algorithms  for  in  situ  location,  by  ultrasound,  of  flaws  in 
plates  or  laminated  layered  materials.  Achieving  this  objec¬ 
tive  requires  detailed  knowledge  of  the  excitation,  propaga¬ 
tion,  scattering  and  detection  of  high  frequency  sound  waves 
in  the  flawed  environment.  Based  on  an  understanding  of  these 
fundamental  wave  phenomena,  one  may  then  attempt  to  construct 
analytical  models  with  accompanying  algorithms,  so  as  to  par¬ 
ametrize  the  NDE  problem  in  terms  of  "good  observables”.  Dur¬ 
ing  the  reporting  period,  carefully  selected  prototype  problems 
have  been  investigated  to  determine  "good  observables"  for  par¬ 
ticular  flawed  environments. 

II .  Accomplishments 

A.  Gaussian  Beam  Excitation  of  a  Solid  Plate 
Near  completion  is  a  study  directed  toward  clarifying  wave 
phenomena  generated  by  a  two-dimensional  Gaussian  beam  input  in 
a  layered  elastic  solid  as  exemplified  (though  the  general  multi¬ 
layer  theory  has  been  developed)  by  an  aluminum  plate  in  vacuum 
[1] .  The  Gaussian  input  has  been  chosen  because  a)  Gaussians  are 
good  models  for  smoothly  tapered  outputs  from  actual  transducers, 
and  b)  they  can  be  used  as  discrete  basis  elements  on  a  (config- 
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uration) -  (spatial  wavenumber)  phase  space  lattice  to  model  rig¬ 
orously  any  transducer  output  [2]  .  The  Gaussian  input  is  anal¬ 
yzed  via  the  complex  source  point  (CSP)  method,  by  which  beam 
excitation  solutions  in  any  (here,  the  layered)  environment  are 
obtained  from  line  or  point  source  excitation  solutions  in  that 
environment  by  assigning  complex  values  to  the  source  coordinate 
location  [3] .  This  avoids  the  need  for  spatial  spectral  decom¬ 
position  of  the  beam  field,  as  is  customary.  Previous  applica¬ 
tions  of  the  CSP  method,  which  requires  careful  study  of  the  an¬ 
alytic  continuation  to  the  complex  source  domain,  have  not  dealt 
with  elastic  layers.  Some  difficulties  encountered  in  this  pro¬ 
cess  are  still  under  investigation,  but  in  the  present  applica¬ 
tion  to  a  compressional  (P)  beam  input,  they  were  de-emphasized 
and  do  not  affect  the  validity  of  the  results.  Extensive  numer¬ 
ical  computations  have  been  performed  to  chart  the  evolution  of 
the  P-beam  input  as  it  interacts  with  the  plate  boundaries  and 
undergoes  P-SV  coupling  at  each  reflection.  The  resulting  dis¬ 
placement  profiles  reveal  clearly  the  beam-like  character  near 
the  source,  the  deterioration  of  the  successively  reflected  beam 
due  to  P-SV  coupling  at  the  boundaries,  and  the  eventual  evolu¬ 
tion  of  oscillatory  mode- like  patterns.  The  results  strongly 
suggest  that  a  hybrid  beam-mode  format  is  worthy  of  considera¬ 
tion  because  the  physical  observables  appear  beam-like  in  cer¬ 
tain  parametric  domains  and  mode- like  in  others. 

B.  Gaussian  Beam  Excitation  of  a  Flawed  Solid  Plate 

In  conjunction  with  the  investigation  of  the  unflawed  plate 
(Section  IIA)  has  been  a  study  aimed  at  determining  "good  para- 
metrizations"  of  observables  when  two-dimensional  Gaussian  P- 
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beam  excitation  of  a  perfectly  bounded  layered  elastic  plate 
(which  appears  like  a  single-layer  thick  plate)  encounters  a 
transversely  thin  region  with  a  weak  bond  imperfection.  The 
flaw  is  assumed  to  give  rise  to  reduced  shear  resistance.  The 
problem  has  been  treated  by  considering  first  the  fields  excited 
by  the  input  beam  in  the  unflawed  environment,  and  then  allowing 
these  fields  to  interact  with  the  flaw;  the  scattered  fields  pro¬ 
duced  thereby  are  monitored  at  the  receiver  on  the  surface.  As 
before,  the  input  beam  was  modeled  compactly  by  a  line  source  at 
a  complex  coordinate  location  (complex  source  point  (CSP)  method 
[3]),  which  permits  the  beam  input  response  to  be  computed  from 
the  line  source  input  response  by  analytic  continuation  of  the 
source  coordinates  into  the  complex  domain. 

The  results  from  the  unflawed  problem  (Section  IIA)  [1]  pro¬ 
vided  the  input  for  the  flawed  configuration.  Since  the  adhe¬ 
sive  layer  was  assumed  to  be  thin  in  comparison  to  the  compres- 

sinal  and  shear  wavelengths,  the  bond  can  be  represented  by  dis- 

*• 

tributed  shear  springs.  Limiting  the  model  to  weak  imperfec¬ 
tions  permits  the  equivalent  surface  fields  induced  on  the  flaw 
to  be  approximated  by  the  incident  field  as  modified  by  the 
stiffness  profile  (Born  Approximation) .  The  scattering  produced 
in  the  plate  by  the  resulting  equivalent  source  distribution  was 
evaluated  by  surface  integration  over  the  flaw  contour,  using 
the  line  force  plate  Green's  function  as  a  kernel.  Primary  at¬ 
tention  was  placed  on  determining  the  physical  observables  (the 
horizontal  and  vertical  displacements)  at  the  detector  on  the  plate 
boundary.  The  parametrization  seeks  to  connect  these  observables 
with  presence,  location,  and  ohysical  constitution  of  the  flawed 
region. 
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In  the  studies  so  far,  the  incident  beam  interacts  with  the 
flaw  before  being  diffused  by  multiple  reflections.  This  makes 
the  beam  perturbed  by  the  flaw  a  "good  observable”  and  suggests  a 
beam  oriented  parametrization.  In  this  endeavor,  a  choice  of 
input  conditions  that  combines  beam  resolution  with  strong  shear 
excitation  at  the  site  of  the  flaw  plays  a  critical  role.  Numer¬ 
ical  results  show  that  a  good  parametrization  can  be  achieved  in 
this  way.  Some  indications  have  been  given  on  how  to  make  the 
inverse  connection  between  surface  displacements  and  flaw  iden¬ 
tification,  but  work  on  this  aspect  is  still  in  progress. 
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INTRODUCTION 

Identification  of  flaws  in  solid  or  laminated  layered  materials  by  ultrasound  requires 
detailed  knowledge  of  the  excitation,  propagation,  scattering  and  detection  of  high 
frequency  sound  waves  in  the  flawed  environment.  An  overall  approach  to  this  problem, 
whereby  one  seeks  to  parametrize  all  relevant  wave  phenomenon  in  terms  of  "good 
observables",  will  be  discussed  elsewhere.  Here,  we  report  on  the  first  phase  dealing  with 
excitation.  In  modeling  actual  or  induced  source  distributions,  we  assign  a  special  role  to 
Gaussian  inputs  because  a)  Gaussians  individually  can  simulate  certain  smoothly  tapered 
distributions,  and  b)  a  superposition  of  Gaussians  on  a  self-consistent  (configuration  - 
wavenumber)  space  lattice  can  represent  any  prescribed  distribution.  To  generate  the 
Gaussian  beams  that  propagate  the  Gaussian  input  into  the  environment,  we  utilize  the 
complex  source  point  (CSP)  method  because  it  avoids  the  plane  wave  spectral  analysis  and 
synthesis  required  in  the  conventional  approach.  By  the  CSP  method.  Green’s  functions 
for  beam  sources  are  derived  directly  from  Green’s  functions  for  line  or  point  sources  by 
analytic  continuation  to  complex  source  locations.  While  the  CSP  method  has  been 
applied  previously  to  beam  propagation  and  diffraction  in  fluid  media  [1-7],  and  to 
studying  such  phenomena  in  semi-infinite  solids  [8-9],  our  concern  here  with  layered 
elastic  media  necessitates  attacking  new  problems  in  the  rigorous  analytic  continuation 
process  that  have  not  been  encountered  in  earlier  studies.  These  problems  are  highlighted 
here,  but  their  resolution  is  still  in  progress  and  will  be  reported  separately.  In  the 
results  cited  below,  excitation  conditions  have  been  chosen  so  as  to  de-emphasize  the 
sources  of  difficulty. 

The  problem  posed  in  the  present  investigation  is  the  mechanism  of  conversion  of  a 
well-collimated  input  beam  in  a  layered  elastic  medium  into  the  spectrum  of  guided 
modes,  favored  among  which  is  a  particular  mode,  one  of  whose  congruences  is  aligned 
with  the  incidence  beam  angle.  The  transitional  phase  is  of  interest  because,  for 
optimum  versatility,  it  may  be  desirable  to  choose  operating  conditions  that  accommodate 
in  a  continuous  manner  beam-type,  mode-type,  or  combinations  of  input  field 
interactions  with  a  fault  zone.  Although  the  formal  theory  has  been  constructed  for  the 
general  case,  we  have  restricted  our  initial  numerical  implementation  to  a  simple  test 
case:  a  P-wave  beam  source  located  in  an  elastic  plate,  which  is  embedded  in  vacuum. 
Results  have  been  generated  via  the  ray-tracing  and  normal  mode  expansion  algorithms, 
published  previously  [10]  but  generalized  here  by  the  CSP  substitution.  Some  of  the 
difficulties  encountered  in  this  generalization  are  summarized  in  the  following  section 
and  typical  numerical  results  are  shown  and  discussed.  The  data  clearly  reveals  the 
influence  of  P-SV  boundary  coupling  on  the  beam-to-mode  conversion  process. 

COMPLEX-SOURCE  POINT  FORMULATION 

To  convert  the  wave  fields  due  to  the  conventional  P-wave  isotropic  line  source 
(Fig.  la)  into  wave  fields  corresponding  to  excitation  by  a  P-wave  Gaussian  beam  source 
(Fig.  lb),  we  shall  employ  the  complex  source  point  technique  mentioned  earlier.  This 
requires  analytic  continuation  of  all  wavefield  expressions  from  the  real  source 
coordinates  (x\  z')  to  the  complex  source  coordinates 

x'  -  x'  +  ib  cos  a  ,  z'  «  z'  +  i  b  sin  a  (I) 


where  (x\  z')  now  locates  t  center  of  the  beam  waist,  b  is  a  bean  1th  parameter,  and 
a  is  the  angle  between  the  beam  axis  and  the  positive  x-axis  (see  Fi*,.  lb  and  Kefs. 

[1-9]).  While  the  analytic  extension  guarantees  formally  that  the  resulting  P  and  SV 
potentials  satisfy  the  appropriate  wave  equations  and  boundary  conditions,  the  validity  of 
a  particular  wavefield  representation  must  be  verified  in  each  case;  i.e.,  the 

representation  obtained  on  replacing  (x',  z)  by  (x\  z)  must  remain  convergent.  It  turns 
out  that  various  alternative  representations  require  different  continuation  strategies. 

When  the  wavefield  representation  is  in  the  form  of  a  spectral  integral  over  the 
spatial  wavenumber  k  corresponding  to  the  x-coordinate,  the  analytic  extension  can  be 
carried  out  systematically  by  contour  deformation,  if  required.  Because  the  integrand 
depends  analytically  on  x'  and  z',  convergence  can  be  tested  continuously  as  x'  and  z'  are 
made  complex.  If,  for  certain  parameter  ranges,  convetgence  at  j  k  |  -»oo  turns  out  to  be 
violated  on  the  original  contour,  that  contour  can  be  deformed  into  appropriate  regions 
of  the  complex  k  plane  wherein  convergence  is  retained.  When  the  representation  is  in 
the  form  of  a  series  like  the  normal  mode  or  the  hybrid  ray-mode  expansion  [10],  direct 

(x'  z')  substitution  may  cause  convergence  problems  due  to  exponentially  growing  factors 
in  the  series  elements.  Contour  deformation  is  now  not  an  immediately  available  option. 
We  shall  report  elsewhere  on  how  to  proceed  in  this  case.  For  presentation  here,  we  have 
chosen  the  beam  parameters  so  that  these  difficulties  can  be  avoided.  We  have  so  far 
dealt  with  the  analytically  extended  ray  and  normal  mode  representations.  Although  both 
formulations  must  (and  do)  generate  the  same  result,  the  ray  formulation  allows  beam 
tracking  through  an  increasing  number  of  reflections  at  the  boundaries  while  the  mode 
representation  highlights  the  fully  excited  guiding  phenomena. 

In  tracking  the  beam  integrals  [6],  it  is  desirable  to  switch  from  the  global  (x,z) 
coordinates  to  the  local  coordinates  defined  with  respect  to  the  axes  of  the  individual 
beams  (see  Fig.  lb).  This  stabilizes  the  analytic  continuation  by  avoiding  oblique  spectra 
with  respect  to  the  beam  axes  at  beam  propagation  angles  approaching  the  x-direction 
and  de-emphasizing  problems  attributable  to  boundary  reflections  of  evanescent  spectra. 
To  illustrate  these  aspects,  we  consider  first  the  incident  P-beam  potential,  which  behaves 
as  though  it  propagates  in  unbounded  space.  From  Eq.  (2)  in  [6j  and  Eq.  (1),  it  is  given 
in  the  global  (x,z)  system  by 

*  =  -  J-  f  —  exp  -f  i[fc_(z-z'-ib  sin  a)+k(x-x'-ib  cos  a)]\dk  (2) 

4tti  -°o  kp  L  J 

where  *p  =  Jkp2-k2,  kp=w/vp,  vp  is.  the  P-wave  propagation  speed,  and  u  is  the  radian 
frequency.  The  integration  path  avoids  the  branch  point  singularities  at  (+kp)  and  (-kp) 
by  indentations  into  the  fourth  and  second  quadrants,  respectively,  of  the  complex  k- 
plane.  The  condition 

z  -  z'  >  b  cos  a  (3) 

guarantees  convergence  of  the  integrand  at  k  — >  ±oo.  In  fact,  the  integrand  decays  in  the 
intervals  from  ±kp  to  too,  as  can  be  verified  by  observing  that  for  k?*^  and  when  Eq.  (3) 
holds,  the  real  part  of  the  phase  in  Eq.  (2)  has  a  negative  derivative  with  respect  to  k. 
This  establishes  the  existence  of  an  evanescent  spectrum.  The  SV-beam  potential 
generated  by  reflection  of  this  incident  beam  from  the  solid-free  space  interface  is,  in 
the  global  (x,z)  system  [6], 

oo  v  r  1 

\j/  m  — L  f  exp  -1  i[*D(z’+ib  sin  a)+*tz+k(x-x'-ib  cos  o))  }>  dk  (4) 

4jtj  j  -oo  Kp  (.  *  ) 

where  r,p  accounts  for  P-S  coupling  at  the  boundary,  k,  -  \/k,2-k*.  k,  -  w/v„  and  v,  is 
the  S-wave  propagation  speed.  The  condition 
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(z'+z)  >  b  cos  at  (jj 

again  ensures  (lie  convergence  of  the  integrand  at  k  — *  ±oo.  Note,  however,  that  the  real 
part  A  of  the  phase  in  Eq.  (4)  is 


A  =  -  v^2-^2  z'  +  kb  cos  a  when  k„>k>kp  (6) 

If 

dA/dk  |  =  —  =*==  z'  +  b  cos  a  >  0  (7) 

*  v4.2-kpJ 

the  spectral  amplitudes  will  grow  in  the  interval  (kp,  k„],  which  behavior  is 
computationally  inconvenient,  and  physically  unnatural  because  the  source  beam  spectrum 
in  that  interval  is  decaying.  The  anomaly  is  avoided  by  choice  of  local  coordinates  and 
corresponding  spectral  decomposition  tied  to  the  incident  beam  so  that  the  beam  angle  in 
this  local  system  equals  a  =  tt/2;  the  local  coordinate  system  is  updated  for  each 
reflection. 

When  the  complex  source  coordinate  substitution  is  performed  on  the  modal  sum, 
the  nonphysical  anomaly  attributed  to  the  kp<k<k,  domain  is  not  avoided  so  easily. 
Pending  further  investigation,  we  have  for  the  present  simply  ignored  the  modes  with 
modal  eigenvalues  km>kp  beca  .se  their  excitation  strength  by  the  specified  P-wave  source 
is  minimal.  Numerical  comparisons  with  fields  computed  by  the  beam  tracking  algorithm 
have  confirmed  the  validity  of  this  assumption. 


A  NUMERICAL  EXAMPLE 


The  mechanism  of  conversion  of  an  initially  well  collimated  incident  beam  in  an 
elastic  plate  into  a  spectrum  of  modes  centered  around  the  one  "best  matched"  to  the 
incident  spectrum  is  illustrated  by  a  numerical  example.  The  parameters  used  in  the 
calculations  pertain  to  an  (aluminum)  plate  (see  Fig.  la): 

vp  =  2.36  x  I06  in./sec.,  v,  =  1.209  x  105  in./sec., 

p  «  2.53  x  10"*  lb.-sec.2/in.4,  a  »  0.1  in.  (8) 

The  60  MHz  P-source  beam  waist  has  been  placed  at  z'  =  0.04  in.,  x'  =  0.0  in.,  with 
beam  width  parameter  b  =  0.1  in.  The  beam  axis  inclination  a  =  0.548  rad.  (see  Fig.  lb) 
lines  up  with  the  upward  P-plane  wave  congruence  of  one  of  the  guided  modes.  The 
axes  of  the  incident  P-beam,  the  reflected  P-bcam  and  the  converted  SV-beams  are 
indicated  in  the  figure.  Fig.  2  shows  the  magnitude  of  the  resulting  P  and  SV  potentials 
$  (dashed)  and  'P  (solid),  respectively,  plotted  in  cross  sections  at  the  axial  distances  x  = 
0.05,  0.08,  0.11,  0.14,  0.17,  0.23  inches  (see  Fig.  lb).  The  potentials  have  been  calculated 
by  the  normal  mode  expansion  in  [10]  and  by  the  ray  series  in  [10].  After  employing  the 
complex  source  coordinates  in  Eq.  (1),  both  procedures  yield  the  same  numerical  values. 

In  the  normal  mode  expansion,  as  mentioned  before,  we  have  retained  only  those  modes 
(64  in  this  example)  with  km<kp  because  physical  considerations  dictate  that  the  incident 
P  beam  couples  very  weakly  to  the  evanescent  P  spectra  contained  in  modes  with  km>kp. 


The  beam-to-mode  conversion  process  is  well  illustrated  by  the  results.  In  the  cross 
section  at  x  »  0.05  in.,  the  incident  P  beam  is  still  far  enough  from  the  upper  boundary 
to  render  the  effect  of  that  boundary  small.  Accordingly,  the  P  potential  retains  its 
essentially  Gaussian  shape  (note  that  the  x  =  constant  planes  do  not  conform  with  the 
natural  beam-centered  coordinate  planes  perpendicular  to  the  beam  axes).  The  ott-axis 
spectra  in  the  Gaussian  tail,  which  do  interact  with  the  boundary,  produce  the  ^ 
decaying  P-SV  reflection  and  the  P-P  reflection  ripple  near  the  top  of  the  P  data.  In  the 
next  observation  plane  at  0.08  in.,  the  incident  beam  axis  has  approached  the  upper 
boundary  sufficiently  closely  to  cause  more  substantial  P-P  reflection  and  P-hv 
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conversion  of  the  olf-axis  ctra.  Inierlcroiico  between  the  uicidc  nd  reflected  P 
spectra  is  evident  by  tlie  more  substantial  ripple,  and  the  converted  .  spectra,  though 
weak  at  the  bottom,  extend  over  the  entire  cross  section.  The  incident  beam  axis  strikes 
the  top  boundary  before  the  next  observation  plane  at  x  =  0.11  in.,  which  therefore 
contains  the  traces  of  both  the  reflected  1*  and  SV  beam  axes.  Strong  interference 
between  the  P  and  P-P  process  near  the  top  boundary  is  evident,  but  one  may  also 
discern  a  P  potential  rise  near  the  bottom,  which  is  due  to  P-SV-P  conversion.  Toe  SV 
spectra  near  the  bottom  boundary  are  now  strong  enough  to  cause  such  conversion,  as 
well  as  the  SV-SV  reflection  ripple  in  the  solid  curve,  whereas  the  bulk  of  the  Sv’ 
reflected  beam  has  moved  from  the  top  toward  the  plate  center.  It  mav  'oe  noted  that  the 
maximum  of  the  wide  SV  reflected  beam  is  shifted  away  from  the  SV  beam  axis 
intercept  with  the  observation  plane.  This  shift  is  due  in  part  to  apparent  asymmetries 
introduced  by  the  substantial  deviation  of  the  x  =  const,  plane  from  the  "natural" 
coordinate  plane  transverse  to  the  beam  axis,  and  also  to  the  lateral  shift  of  the  reflected 
SV  beam  axis  on  the  top  boundary  which  is  known  to  accompany  the  P-SV  conversion. 
These  phenomena  are  further  enhanced  at  x  =  0.14  in.,  and  also  at  x  =  0.17  in.  where 
substantial  reflection  and  P  conversion  of  the  downgoing  SV  beam  at  the  bottom  make 
their  contribution.  Clearly,  individual  beams  (certainly  for  SV  waves  but  also  for  P 
waves)  are  now  difficult  to  resolve,  and  the  potentials  begin  to  stabilize  toward 
oscillatory  patterns  with  fixed  periodicity.  These  periodicities  turn  out  to  be  essentially 
those  found  in  the  guided  mode  whose  upgoing  plane  wave  P  spectra  line  up  most  closely 
(in  our  case,  exactly)  with  the  dominant  wave  spectra  of  the  incident  P  beam  (i.e.,  with 
the  beam  axis  if  the  incident  beam  is  well  collimated).  The  SV  periodicity  pattern  is 
already  well  established  at  x  =  0.17  in.  because  the  SV  beams  with  their  steeper 
propagation  angles  have  had  an  opportunity  to  interact  with  both  the  top  and  bottom 
surfaces  of  the  plate.  For  the  P  patterns,  larger  x  values  are  required  before  the 
analogous  interactions  and  consequent  periodicity  stabilizations  occur;  we  have  verified 
that  x  =  0.23  in.  is  adequate  (see  Fig.  2).  Other  more  extensive  calculations  are  in 
progress  and  will  be  reported  elsewhere. 

DISCUSSION 

In  this  paper,  we  have  explored  the  mechanism  by  which  an  incident  well 
collimated  beam  excites  the  spectrum  of  guided  modes  that  can  be  supported  in  a  layered 
structure.  For  simplicity,  we  have  chosen  a  single  elastic  plate  embedded  in  vacuum,  but 
other  models,  with  exterior  fluid  loading,  are  already  under  consideration.  Results  for 
the  P  and  SV  potentials  in  successive  cross  sections  have  been  computed  numerically  by 
the  beam  tracking  algorithm  and  by  normal  mode  expansion.  Both  being  formally  exact, 
they  are  expected  to  (and  they  do)  furnish  the  same  numerical  data.  Concerning  the 
form  of  the  result,  the  field  in  initial  cross  sections  is  beam  shaped,  reproducing 
essentially  the  incident  beam  in  unbounded  space,  before  it  reaches  one  of  the 
boundaries.  This  suggests  that  the  wave  fields  in  the  plate  should  be  calculated  in  this 
parameter  regime  by  beam  functions.  Expansion  in  terms  of  the  spectrum  of  guided 
modes  in  the  plate  is  equally  valid,  but  many  modes  are  required  to  synthesize  the 
Gaussian  beam  profile.  Thus,  guided  modes  are  phenomenologically  unsuited  and  furnish 
poor  observables  here,  whereas  beams  furnish  good  observables.  After  the  beam  has 
interacted  with  one  of  the  boundaries,  the  resulting  interference  and  conversion 
phenomena  blur  the  originally  clean  beam  profile,  and  after  successive  reflections,  any 
resemblance  to  beams  has  disappeared  altogether.  The  wave  pattern  in  the  cross  section 
becomes  oscillatory,  thereby  implying  that  the  good  observables  are  now  the  guided 
modes.  This  suggests  use  of  our  beam-generated  self-consistent  hybrid  ray-mode 
algorithm  [10],  which  can  be  structured  so  that  it  encompasses  automatically  the  transition 
from  one  good  observable  parametrization  to  the  other.  This  aspect  is  now  under 
investigation. 

One  interesting  observation  lias  already  emerged  from  these  numerical  results.  It 
was  noted  earlier  that  a  few  reflections  between  the  plate  boundaries  suffice  to  diffuse 
the  profile  of  an  initially  well  collimated  beam  into  a  cross  sectional  standing  wave 
pattern.  The  periodicity  of  this  pattern  is  well  reproduced  by  that  of  the  guided  mode, 
one  of  whose  constituent  plane  wave  congruences  is  closely  aligned  with  the  dominant 
wave  spectra  that  propagate  along  the  axis  of  the  incident  beam.  Thus,  the  total  field 
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stabilizes  rather  quickly  to  j  the  periodicity  pattern  of  the  mode  ored  by  the 
incident  beam.  However,  that  mode  is  not  excited  in  pure  form  as  v.  .denced  by  the 
nonuniform  amplitude  profile  of  the  total  field.  Therefore,  the  function  of  modes  with 
other  periodicities  seems  to  be  primarily  to  adjust  the  amplitude:  in  the  oscillatory 
pattern  to  their  proper  values.  The  close  matching  of  the  total  field  periodicities  with 
those  of  the  favored  mode  is  important  because  good  prediction  of  null  locations  is 
required  for  wavefield  interaction  with  structural  faults.  Further  study  of  this 
phenomenon  is  intended  to  ascertain  whether  these  observations  apply  also  under  more 
general  conditions. 
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(a)  real  source  coordinates  (x\  z')  (P-wave  isotropic  line  source). 
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(b)  complex  source  coordinates  (x\  ?)  (P-wave  beam  source).  The  beam  parameter  b 
in  (1)  is  related  to  the  (l/e)  beam  width  w,  at  the  waist,  we  =  (2b/k  )*/J.  Also 
shown  on  this  figure  are  the  axes  of  the  incident  and  reflected  P-beams  (dashed 
lines),  and  of  the  coupled  SV- beams  (solid  lines),  and  their  relationship  to  the 
cross  sections  in  Fig.  2. 


Fig.  1.  Elastic  plate  with  thickness  a,  characterized  by  wave  velocities  vp,  vt  and  by 
density  p. 


Fig.  2.  Amplitude  profiles  of  the  P-wave  potential  $  (dashed)  and  the  SV-wave 

potential  (solid)  due  to  a  P-beam  input,  observed  at  successive  plate  cross 
sections  in  the  interval  from  x  =  0.05  in.  to  x  =  0.23  in.  The  plate  parameters 
are  in  (8)  and  the  frequency  is  60  MHz.  The  beam  parameters  are:  x'  »  0,  z'  = 
0.04  in.,  a  -  0.548,  b  =  0.1  in.  The  vertical  axis  measures  the  cross  sectional 
coordinate  z.  The  horizontal  axis  measures  |  $  |  and  |  9  | ,  respectively,  on 
different  scales,  with  the  maximum  in  each  wave  denoted  by  the  number 
listed.  This  figure  should  be  viewed  together  with  Fig.  lb,  which  shows  the 
disposition  of  the  incident  and  reflected  beam  axis. 


Appendix  B 

To  be  pubx-dhed  in  the  Proceedings  of  the  ^  .AM-IUPAP 
Symposium  on  Elastic  Wave  Propagation,  held  in  Galway, 
Ireland  20-25  March,  1988 
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Location  and  Identification  of  faults  In  multilayer  elastic  materials  by 
ultrasound  Is  aided  by  a  physically  based  parametrlzatlon  of  the  Input, 
scattered  and  detected  fields.  When  the  transducer  Input  is  beam-shaped, 
the  beam-to-mode  conversion  in  the  unflawed  layered  environment  sug¬ 
gests  a  "good"  parametrization  in  terms  of  a  self-consistent  hybrid  beam¬ 
mode  format.  The  scattered  field  produced  by  interaction  of  this  beam¬ 
mode  field  with  a  fault  zone  should  then  be  parametrized  In  a  similar 
manner.  This  strategy  guides  the  present  Investigation  of  a  weak  bonding 
flaw  in  a  multilayer  aluminum  plate.  The  horizontal  and  vertical  displace¬ 
ments  excited  by  a  high  frequency  two-dimensional  dllatatlonal  (P)  Gaus¬ 
sian  input  beam  have  previously  been  tracked  through  successive  cross 
sections  in  the  perfectly  bonded  material.  The  resulting  displacement 
profiles  reveal  clearly  the  beam-Uke  character  near  the  source,  the 
deterioration  of  the  successively  reflected  beam  due  to  P-SV  coupling  at 
the  boundaries,  and  the  eventual  evolution  of  oscillatory  mode-like  pat¬ 
terns.  This  Input  is  now  allowed  to  interact  with  an  elongated  weak  bond 
zone.  The  induced  equivalent  forcing  terms  are  modeled  In  the  Born 
approximation,  and  the  scattered  field  Is  evaluated  accordingly.  Depend¬ 
ing  on  the  flaw  size,  its  location  relative  to  the  Input  and  output  transduc¬ 
ers,  and  other  variables,  the  detected  response  at  the  plate  surface  may 
contain  beam-like  or  mode-like  features.  The  beam-Uke  phenomena  are 
explored  here  with  a  view  toward  finding  conditions  through  which  the 
physical  observables  that  should  facilitate  flaw  location  and  identification 
are  enhanced.  Although,  for  convenience,  the  numerical  data  have  been 
generated  by  normal  mode  summation,  the  results  reveal  clearly  that  the 
hybrid  beam-mode  format,  to  be  developed  next,  furnishes  the  proper 
parametrization. 


1.  INTRODUCTION 


Analytical  modeling  of  flaw  identification  in  a  layered  composite  through  response  to 
an  ultrasonic  beam  input  poses  a  problem  of  substantial  complexity.  A  good 
identification  scheme  must  rely  on  a  good  parametrization  of  the  propagation  and 
scattering  process  in  terms  of  physical  observables.  Translating  such  observables  into 
analytic  form  is  usually  done  best  by  selecting  coordinates  that  are  well  matched  to  the 
phenomenology.  This  is  easily  accomplished  when  considering  the  input  beam,  the 
flaw  induced  scattering,  and  the  layered  environment  in  Isolation  but  it  is  far  less  tran¬ 
sparent  when  all  of  these  constituents  interact:  the  beam  coordinates  are  generally  not 
matched  to  those ^pf  the  flaw  nor  to  thqpe  of  the  layers.  Thus,  analytical  basis  eie- 


ments  that  describe  compactly  the  wave  phenomena  associate  with  one  of  the  consti¬ 
tuents  are  likely  to  be  awkward  and  less  physical  for  the  other  constituents.  One 
attempt  to  cope  with  these  difficulties  is  via  self-consistent  hybrid  combinations  of 
different  basis  elements,  chosen  so  as  to  mode!  the  various  observables  in  those 
parameter  ranges  where  they  are  most  pronounced. 


In  a  previous  and  the  present  Investigation,  an  effort  has  been  made  to  address  these 
concerns  systematically.  The  overall  problem  has  been  decomposed  by  considering 
llrsl  the  fields  excited  by  the  Input  beam  In  the  unflawed  environment,  to  be  followed 
by  the  interaction  of  these  fields  with  the  flaw,  and  the  consequent  scattered  fields 
produced  thereby  at  the  receiver.  For  simplicity,  the  problem  is  taken  to  be  two- 
dimensional,  and  the  unflawed  material  is  a  perfectly  bonded  layered  aluminum  plate 
(which  therefore  acts  like  a  single  thick  layer)  In  vacuum.  The  input  beam  Is  assumed 
to  be  compresslonal  (P-wave)  and  to  have  a  Gaussian  profile  so  that  it  can  be 
modeled  compactly  by  a  line  source  at  a  complex  coordinate  location  (complex  source 
point  (CSP)  method)  [ l] .  This  permits  the  beam  input  response  to  be  computed  from 
the  line  force  input  response  by  analytic  continuation  of  the  source  coordinates  into 
the  complex  domain.  Depending  on  whether  the  observer  is  near,  or  far  from,  the 
source  region,  the  field  In  the  plate  has  a  beam-like  or  mode-like  character,  undergo¬ 
ing  rapid  conversion  from  the  former  to  the  latter  due  to  multiple  reflections,  with 
consequent  compressional-shear  (P-SV)  coupling,  at  the  plate  boundaries.  A  hybrid 
ray-mode  form  can  describe  this  phenomenology.  This  first  (unflawed)  phase  of  the 
problem  has  already  been  treated  In  detail  [2],  and  the  results  from  that  Investigation 
provide  the  Input  for  the  flawed  environment  here. 


To  test  the  proposed  problem  strategy  without  undue  analytic  complication,  an  imper¬ 
fect  bond  or  finite  extent  aligned  parallel  to  the  plate  boundaries  has  been  postulated. 
The  imperfection  pertains  only  to  reduced  shear  resistance,  and  thus,  since  the 
adhesive  layer  is  thin  In  comparison  to  the  wavelengths,  the  bond  can  be  represented 
by  distributed  shear  springs.  For  weak  imperfections,  the  equivalent  surface  fields 
induced  on  the  flaw  can  be  approximated  by  the  Incident  field  as  modified  by  the 
stiffness  profile  (Born  approximation);  because  the  flaw  responds  to  tangential  shear 
but  not  to  compression,  the  Incident  field  at  the  flaw  site  should  have  a  strong  longitu¬ 
dinal  shear  component.  The  scattering  produced  in  the  plate  by  the  resulting 
equivalent  source  distribution  "can  be  evaluated  by  surface  integration  over  the  flaw 
contour,  using  the  line  force  plate  Green's  function.  The  detector  Is  assumed  to  be 
located  on  the  plate  boundary;  therefore,  primary  attention  is  placed  on  determining 
the  physical  observables  (i.e.,  the  horizontal  and  vertical  displacements)  there. 


The  presentation  below  begins  with  a  summary  description  of  beam  excited  fields  that 
have  favorable  properties  for  Interaction  with  the  flaw  in  the  thick  plate,  and  of  good 
parametrizations  of  these  fields.  The  surface  integral  for  the  scattered  field  Is  formu¬ 
lated  next,  with  expressions  given  for  the  scattering  induced  displacements  on  the 
plate  surface.  The  detection  strategy  emphasizes  beam-like  phenomena.  Qualitative 
measures  of  "good  observables"  are  Introduced,  which  permit  drawing  inferences 
about  the  flaw  size  and  location  from  the  detected  data.  Bach  of  these  phases  is 
accompanied  by  numerical  calculations  that  demonstrate  the  detailed  field  behavior. 
The  paper  ends  with  conclusions  based  on  these  results. 


2.  FORMULATION 
2.1  General  Concepts 


The  physical  configuration,  schematized  in  Fig.  1,  involves  a  bonded  plate  of  thickness 
a,  constitutive  parameters  (Lame*  constants)  X,  n  and  (density)  p,  with  excitation  pro- 


vided  by  a  two-dimens  iv/  nal  ultrasonic  harmonic  Gaussian  P-beam  input  whose  waist  is 
centered  at  (x,z)  =  (0,z')  and  whose  inclination  with  respect  to  the  x-axis  is  a.  A 
longitudinally  oriented  thin  boundary  flaw  of  length  2# Is  centered  at  (xt,  zf).  As  noted 
in  Section  l,  the  bond  imperfection  Is  only  with  respect  to  its  shear  resistance,  while 
its  full  strength  is  retained  in  the  normal  direction.  The  rigidity  of  the  weakened  bond 
is  represented  by  a  distribution  of  equivalent  springs  along  the  flaw  contour,  the 
equivalent  spring  stiffness  being  defined  as 

K=rzx/Au  (1) 

where  rzx  is  the  shearing  stress  at  the  Joint,  and  Au  =  u+-  u"  denotes  the  x- 
displacement  Jump  discontinuity  across  the  bond  line. 
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Fig.  l.  Gaussian  P-beam  input  to  elastic  plate  with  bonding  flaw  of 

length  2$  centered  at  (x,,zf).  Incident  and  reflected  P-beam 
axes  (dashed)  and  reflected  SV-beam  axis  (solid)  are 
included.  Numerical  data  in  Figs.  2-8  are  for  an  aluminum 
plate  of  thickness  a  =  0.3  in.  and  bond  parameters  x,= 

0.105,  0.135,  0.09,  0.045,  Z{=0.15  in.,  2)2=  0.090  in.  The 
beam  is  modeled  by  the  complex  source  point  method,  with 
complex  source  point  location  x'  =  lb  cos  a,  z'= 
z'  -fib  sin  a,  z' =  0.08  in.  ,  a=45°,  b  =  0.3  in.  The 
corresponding  beam  waist  is  we  =  (2b/kp)1/2  =  0.0335  in. 

The  objective  is  the  calculation  of  the  horizontal  and  vertical  displacements  u(x,z)  and 
w(x,z),  respectively,  which  constitute  the  physical  observables  generated  by  input 
forcing  functions.  Because  the  problem  is  composite  in  the  sense  described  in  Section 
1,  various  analytical  building  blocks  are  employed  to  characterize  the  relevant  wave 
phenomena  in  a  manner  amenable  to  good  Interpretation.  This  may  require  use  of 
different  coordinate  systems  "matched"  to  particular  phenomena  (for  example,  beam 
centered  coordinates  instead  of  the  plate  coordinates  (x,z)).  For  the  force  excited  dis¬ 
placements,  it  may  be  convenient  to  use  the  scalar  potentials  $(x,z)  and  'p(x.z)  per¬ 
taining  to  compresslonal  (P)  and  vertically  polarized  (SV)  motion,  respectively,  with 
propagation  speeds  vp  =  ((X+aO/p]  1/2  and  vs  =  Im/p]  1/a-  The  displacements  are 
derived  from  the  potentials  via  the  differentiations 


d<$  _  d^r  _  d$  d4i 

dx  dz  W  dz  dx 


(2) 


Wave  phenomena  In  the  composite  environment  of  Fig.  1  can  be  synthesized  by  spec¬ 
tral  decomposition  and  recombination  with  respect  to  relevant  "preferred"  coordinates. 
For  each  individual  beam,  whether  Incident  or  multiply  reflected,  this  is  best  done  in 
the  plane  perpendicular  to  the  beam  axis.  For  waves  “tied”  to  the  plate,  the  spectral 
decomposition  is  along  x  or  z.  Decomposition  along  x  organizes  the  response  in 
terms  of  a  plane  wave  spectral  continuum  with  horizontal  wavenumbers  k, 

00 

A(x,z)  =  -i-  /  A(z,/c)  ellcx  dx  (3) 
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where  A(z ,«),  with  k  =  /c(k)  representing  the  vertical  wavenumber,  is  the  z- 
dependent  spectral  wave  amplitude  corresponding  to  the  physical  quantity  A(x,z). 


Decomposition  along  z  organizes  the  response  In  terms  of  a  discrete  spectrum  ol' 
guided  modes  with  vertical  wavenumber  «m, 

.  ikm(Km)x 


A(x,z)  =  £  Am(z,/cm)  e 
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(4) 


where  m  is  the  mode  index  and  Am(z,Km)  is  the  mode  function.  The  wavenumbers 
k  and  k  are  related  by  the  dispersion  equations 

s)2  for  P,  SV  waves,  (5) 


+  *p,s  =  (Wv 


P,S/ 


with  lm(kor/c)>0  to  ensure  convergence  when  the  wavenumber  Is  nonreal.  The 
continuous  wave  functions  A  In  (3)  or  each  mode  function  Am  In  (4)  can  be 
expressed  either  in  oscillatory  form  or  In  terms  of  upgoing  and  downgoing  plane  wave 
congruences  with  dependence  exp(±i/cz)  or  exp(±l«mz).  The  latter,  more  general, 
format  identifies  plane  waves  as  the  fundamental  spectral  building  blocks.  It  should 
be  noted  that  P-SV  coupling  at  the  layer  boundaries  implies  that  both  wave  types  have 
the  same  horizontal  wavenumber  k. 


The  spectral  considerations  for  the  beam-excited  layer  have  been  described  in  detail  in 
a  separate  publication  (2].  The  phenomenological  summary  above  has  been  presented 
to  provide  Justification  for  subsequently  used  qualitative  estimates  that  explain  the 
numerically  generated  displacement  observables.  For  example,  the  displacements  in 
(2)  involve  contributions  from  the  differentiated  compresslonal  and  shear  potentials 
simultaneously,  thereby  mitigating  against  the  possibility  of  unscrambling  their 
separate  impact  in  the  observed  waveforms.  However,  In  either  the  x-based  or  z- 
based  spectral  domains,  it  follows  for  each  spectral  plane  wave  element 
~  exp(i/cz  4-  ikx)  that  the  derivative  operations  imply 

•#--w  .4--*  ■'  <s> 
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Thus,  depending  on  the  ratio  >c/k  (i.e.,  the  spectral  plane  wave  propagation  angle), 
one  may  encounter  parameter  fequlrements  where  either  <I>  or  ^  predominates.  The 
observed  displacements  can  then  be  Interpreted  from  the  form  of  the  appropriate 
dominant  potential.  Conversely,  this  guideline  has  been  employed  by  us  in  the  selec¬ 
tion  of  Input  beams. that  favor  the  excitation  of  strong  compression  or  strong  shear. 


2.2  Implementation 


2.2.1  Beam  excited  layer  without  flaw 


The  solution  strategy  for  charting  the  evolution  of  an  Incident  Gaussian  P-wave  beam 
in  the  perfectly  bonded  layer  has  been  treated  analytically  and  numerically  In  reference 
(2],  The  beam  input  has  been  generated  from  the  line  force  input  by  analytic  con¬ 
tinuation  of  the  source  point  coordinates  from  real  to  complex  values,  and  the  poten¬ 
tials,  stresses  and  displacements  have  been  calculated  In  various  vertical  cross  sections 
and  also  along  the  surface.  Guided  by  the  spectral  concepts  discussed  In  Section  2.1, 
it  has  been  possible  to  explain  the  observed  waveforms  In  terms  of  multiply  reflected 
P  and  SV  beams  near  the  Input,  and  in  terms  of  the  mode  "most  favored"  by  the 
incident  beam  la  the  Intermediate  and  more  distant  regions.  It  has  also  been  possible, 
based  on  these  spectral  considerations,  to  "tune”  the  Input  beam  parameters  so  as  to 
generate  a  desired  response,  for  example,  strong  horizontal  shear  In  the  longitudinal 
plane  that  Is  to  contain  the  flaw.  A  typical  set  of  response  curves  Is  shown  in  Fig.  2. 


2.2.2  Excitation  of  the  flaw 


With  the  flaw  in  place,  as  shown  in  Fig.  l,  and  knowing  the  incident  field  at  the  flaw 
from  the  results  in  Section  2.2.1,  one  may  formulate  the  scattering  problem.  In  gen¬ 
eral,  this  requires  solution  of  a  surface  integral  equation  for  the  unknown  induced 
sources  on  the  flaw.  However,  in  view  of  the  assumed  weakness  of  the  imperfection, 
the  induced  sources  can  be  approximated  by  those  corresponding  to  the  incident  field 
on  the  flaw  (Born  approximation).  Accordingly,  from  (l),  the  displacement  discon¬ 
tinuity  across  the  weakened  zone  becomes 

Tx 

Auso  =  —  (7) 

where  superscripts  "i"  and  ''sc"  refer  to  incident  and  scattered  quantities,  respectively. 


3.  SCATTERED  FIELDS 


The  displacements  usc  and  wsc  generated  by  the  initial  conditions  In  (7)  can  be 
evaluated  by  use  of  Green's  function  techniques  [3], 

usc(x,z)  =  /  £zxx(x',zf;x,z)Ausc(x\z,)  dx' 

‘r-j? 

w3C(x,z)  =  /  Ezxz  (x(zf;x,z)Au*c(x',zf)dx'  (8) 

The  Green's  stress  functions  Ezxx  and  Ezxz  in  (8)  permit  calculation  of  the  shearing 
stress  rzx(x,z()  in  a  perfectly  bonded  layer  due  to  a  horizontally  or  vertically  polarized 
input  force,  respectively.  They*-can  be  represented  in  any  of  the  alternative  forms  dis¬ 
cussed  in  reference  [4,5].  In  the  calculations  that  follow,  the  modal  expansion  form 
was  chosen  for  convenient  reference,  but  the  interpretation  of  the  results  has  been  car¬ 
ried  out  by  recourse  to  spectral  concepts  tied  to  "good  observables,"  as  described  in 
Section  2.1. 


4.  NUMERICAL  RESULTS 


An  extensive  sequence  of  numerical  calculations  has  been  carried  out  to  gain  under¬ 
standing  of  the  fields  generated  in  the  composite  geometry  of  Fig.  1,  in  particular,  the 
displacements  on  the  upper  boundary  of  the  layer,  where  the  detector  is  assumed  to 
be  located.  To  this  end,  the  displacements  generated  by  the  input  beam  without  the 
flaw,  by  the  scattering  in  (8),  and  by  the  combination  of  both,  have  been  calculated 
and  plotted  individually.  An  effort  has  been  made  to  choose  the  input  beam  so  as  to 
give  rise  in  the  total  displacements  to  distinct  features  that  can  be  tied  to  flaw  location 
and  strength.  This  has  led  to  oblique  input  beams  that  generate  resolvable  multiple 
excursions  over  the  extent  of  the  flaw.  The  scattered  shear  waves  then  become 
clearly  visible  on  the  surface  even  in  the  total  response.  Deconvolution  techniques  for 
separate  extraction  of  the  scattered  field  displacements  can  furnish  an  even  more  use¬ 
ful  data  base. 


Figures  2-0  contain  the  numerical  data  which  have  been  selected  for  presentation 


here.  The  parameters 


.id  pertain  to  an  aluminum  plate  (s^e  Fig.  l): 
vp  =  2.36x10s  in./sec.,  vs  =  1.209x10s  in./sec. 

p  =  2.53x10” 4  lb.-sec.2/ln.\  a  =  0.3  in.  (9) 


The  flaw  length  is  2$  =  0.000  in.,  and  the  shear  spring  stiffness  distribution  is  assumed 
to  be  Gaussian 


K  =  K0e[(x  ,  Kc  =  5x  109  ib./in.3  . 


The  bond  location  at  z(  =  0.150ln.  is  midway  between  the  plate  boundaries,  and  the 
center  of  the  flaw  for  the  four  cases  considered  is  located  at  xf  =  0.165,  0.135,  0.090, 
0.045  in. 
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Fig.  2. 

Shearing  stress  generated  along  bond 
line  of  unflawed  plate.  Note  that  the 
first  peak  corresponds  to  the  incident 
P-beam  and  the  second  peak  to  the 
converted  SV-beam  (see  Fig.  1).  A 
number  next  to  a  peak  denotes  the 
peak  amplitude. 
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Fig.  3. 

Induced  source  distribution  for  four 
flaw  locations.  Gaussian  spring  stiffness 
distribution  K  =  K„  exp[(x-  x,)/$  ]2, 
K0  =  5xl09  lb./in3.  Asymmetries  in  the 
distribution  for  Cases  1,  2,  4  arise 
because  the  incident  beam  illumination 
over  the  flaw  is  asymmetrical. 


The  waist  of  the  20MHz  P-source  beam  has  been  placed  at  x'  =  0.0  in.,  z'  —  O.oein. 
with  beam  width  parameter  b  =  0.3  in.  and  corresponding  beam  waist 
we  =  (2b/kp)1/2  =  0.0335  in.  The  beam  axis  Inclination  angle  a  =  45°  accommodates 
the  resolvability  requirement  mentioned  above,  and  It  also  leads  to  excitation  of 
strong  horizontal  shear  at  the  flaw  site.  For  these  Incident  beam  conditions,  there  Is  a 
cluster  of  modes  (with  P-wave  congruences  surrounding  the  45°  axis)  which  contri¬ 
butes  appreciably  to  the  shear  Inducing  field  on  the  flaw;  only  antisymmetric  modes 
need  to  be  included  for  the  calculations  because  symmetric  modes  have  no  horizontal 
shear  at  that  location.  Out  of  the  contributing  spectrum  of  64  modes  (see  [2]),  a  total 
of  20  (Ailing  the  angular  spectrum  Interval  58°<0<3O°)  was  retained  for  the  calcula¬ 
tions;  the  error  produced  by  omitting  the  others  is  too  small  to  be  discerned  on  the 
plots.  The  P-wave  angular  spectrum  spread  Is  compatible  with  the  important  spectral 
content  of  the  incident  beam.  Spectral  filtering  tests  have  been  performed  to  assess 
the  loss  of  definition  when  the  mode  spectrum  is  truncated  further.  For  example, 
retaining  only  7  modes  in  the  cluster,  which  fill  the  angular  spectrum  interval 
38°<0 <52°,  one  finds  a  maximum  amplitude  error  of  1%  but  more  significant  error 
in  the  phase.  By  choosing  a  larger  beam  incidence  angle,  one  may  excite  a  particular 
mode  in  purer  form  because  the  mode  separation  increases  under  these  conditions. 
However,  this  is  accompanied  by  a  reduction  of  the  inducing  shear,  and  also  by  poorer 


rcsolvabiliiy  of  the  mu"  ->le  reflected  beams.  Such  a  strate  may  be  useful  for  longi¬ 
tudinally  extended  bo,  with  strong  imperfections  which  .an  appreciably  aiTect  the 
mode  structure  in  the  hawed  interval,  thereby  making  the  perturbed  mode  a  good 
observable. 


Fig.  4.  Scattered  vertical  and  horizontal  displacement  fields  on  upper 

surface  of  plate.  First  peaks  correspond  to  intersection  of 
direct  scattered  SV-beam  with  the  surface. 


All  numerical  calculations  were  performed  on  a  PC/AT  compatible  Everex  10MHz 
System  1800,  having  a  FAST  10MHz  80287  Math  Compressor.  The  scattered  dis¬ 
placement  fields  were  determined  by  evaluating  the  Integrals  In  (8)  numerically,  using 
30  nodal  points  (5  to  8  nodes  per  wavelength)  along  the  flaw.  The  shearing  stress 
generated  by  the  Input  beam  along  the  bond  line  of  the  unflawed  plate  Is  shown  in 
Fig.  2.  The  axes  of  the  Incident  beam,  the  reflected  P-beam  and  the  converted  SV- 
beam  indicate  that  the  first  peak  corresponds  to  the  Incident  P-beam,  while  the  second 
peak  Is  associated  with  the  converted  SV-beam.  For  a  steeper  beam  angle  a,  the  hor¬ 
izontal  shearing  stress  due  to  the  P-beam  is  decreased,  rendering  the  reflected  SV  con¬ 
tribution  dominant.  The  Induced  sources  along  the  flaw,  computed  from  (7),  are 
given  in  Fig.  3  for  the  four  flapr  locations.  For  Case  3,  where  the  beam  axis  strikes 
the  center  of  the  flaw,  the  Induced  source  magnitudes  are  largest,  and  Gaussian  sym¬ 
metry  is  closely  maintained.  The  half-power  incident  beam  width  projected  onto  the 
flaw  site  is  0.05  in. 


The  scattered  displacement  fields  on  the  upper  surface  of  the  plate  are  shown  in  Figs. 
4  for  the  four  flaw  locations.  The  scattered  shear  waves  become  clearly  visible,  with 
details  for  Case  3  gleaned  from  Fig.  5.  Note  that  the  first  peak  occurs  at  the  intersec¬ 
tion  of  the  direct  scattered  SV-beam  and  the  surface,  and  that  the  shape  about  that 
peak  has  similarities  with  the  shape  of  the  induced  source  distribution.  Subsequent 
peaks  line  up  with  the  multiple  reflected  SV-beams. 


The  displacements  on  the  surface  generated  by  the  input  beam  In  the  plate,  by  the 
scattering  from  the  flaw,  and  by  the  combination  of  both  have  been  plotted  Individu¬ 
ally  In  Fig.  6  for  Case  3.  Clearly  visible  In  a  region  which  was  Initially  quiescent  is  the 
disturbance  attributable  to  the  scatterer.  The  corresponding  real  displacement  com¬ 
ponents  are  also  shown.  The  observed  periodicities  turn  out  to  be  essentially  those 
found  In  the  guided  mode  whose  upgoing  reflected  plane  wave  SV-spectra  line  up  best 
with  the  dominant  wave  spectra  of  the  incident  P-beam;  the  amplitude  modulation  is 
due  to  interference  between  several  of  the  modes  in  the  surrounding  cluster.  The 
small  oscillatory  precursor  in  the  total  displacements,  which  appears  to  the  left  of  the 
undisturbed  region  in  the  unflawed  plate  response,  Is  explicitly  tied  to  the  flawed 
environment. 
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ginary  components  (not  shown)  exhibit 
similar  characteristics.  Phase  plots  (not 
shown)  possess  essentially  the  linear 
variation  associated  with  the  P-exclted 
incident  SV-plane  wave  observed  at  the 
boundary,  except  near  the  onset  of  the 
first  disturbance  where  the  field  departs 
strongly  from  the  plane  wave  prototype. 

To  link  the  measured  displacements  on  the  surface  with  the  flaw  that  induces  the  per¬ 
turbations,  we  are  continuing  the  search  for  other  observables,  in  addition  to  those 
described  above.  There  are  Indications  that  the  readily  accessible  ratio  (u/w)  contains 
information  which  can  be  related  to  the  flaw  location  by  simple  ray  construction,  while 
inferences  on  flaw  strength  can  be  drawn  from  the  data  for  u  or  w  individually. 
Details  of  the  inversion  process,  i.e.,  specifying  flaw  parameters  from  measured  dis¬ 
placements  at  the  surface,  are  still  under  Investigation  and  will  be  reported  separately. 
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5.  SUMMARY 


This  study  has  been  directed  toward  a  "good  parametrlzatlon"  of  observables  when  a 
Gaussian  P-beam  input  in  a  perfectly  bonded  layered  two-dimensional  elastic  plate 
encounters  a  transversely  thin  region  with  weak  debonding  of  shearing  stress.  The 


observables  are  the  horizontal  and  vertical  displacements  on  the  surface  of  the  plate; 
the  parametrization  seeks  to  connect  these  observables  with  presence,  location,  and 
physical  constitution  of  the  flawed  region.  Emphasis  here  has  been  on  how  the  plate 
environment,  first  without  and  then  with  the  flaw,  affects  the  incident  beam;  i.e.,  the 
parametrization  is  beam  oriented.  In  this  endeavor,  the  choice  of  input  conditions 
that  combines  multiple  beam  resolution  with  strong  shear  excitation  at  the  site  of  the 
flaw  plays  a  critical  role.  The  numerical  results  have  shown  that  a  good  parametriza¬ 
tion  can  be  achieved.  Some  Indications  have  been  given  on  how  to  make  the  inverse 
connection  between  surface  displacements  and  flaw  Identification  but  work  on  this 
aspect  is  still  in  progress. 


Despite  the  beam  oriented  parametrization,  the  results  have  been  computed  by  mode 
summation  because  the  complex-source-point  adapted  mode  algorithm  has  been  avail¬ 
able  from  a  prior  phase  in  this  investigation  [2].  Evidently,  the  numerical  algorithm 
should  be  restructured  in  hybrid  beam-mode  form  (2,5,6)  to  highlight  the 
parametrized  observables  directly.  Consideration  of  this  aspect  is  postponed  until  a 
thorough  understanding  has  been  gained  of  the  variety  of  phenomena  that  accompany 
the  flaw  identification  problem  in  a  multilayer  bonded  plate  environment. 
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